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An interacting one-dimensional electron system, the Luttinger liquid, is distinct from the "con- 
ventional" Fermi liquids formed by interacting electrons in two and three dimensions Some of 
its most spectacular properties are revealed in the process of electron tunneling: as a function of the 
applied bias or temperature the tunneling current demonstrates a non-trivial power-law suppression 
@, l|- Here, we create a system which emulates tunneling in a Luttinger liquid, by controlling the 
interaction of the tunneling electron with its environment. We further replace a single tunneling 
barrier with a double-barrier resonant level structure and investigate resonant tunneling between 
Luttinger liquids. For the first time, we observe perfect transparency of the resonant level embedded 
in the interacting environment, while the width of the resonance tends to zero. We argue that this 
unique behavior results from many-body physics of interacting electrons and signals the presence of 
a quantum phase transition (QPT) 4J. In our samples many parameters, including the interaction 
strength, can be precisely controlled; thus, we have created an attractive model system for studying 
quantum critical phenomena in general. Our work therefore has broadly reaching implications for 
understanding QPTs in more complex systems, such as cold atoms [f| and strongly correlated bulk 
materials [f|. 



Unlike two- and three-dimensional Fermi liquids, a 
Luttinger liquid completely "dissolves" individual elec- 
trons, replacing them with collective plasmon waves. 
When in the process of quantum-mechanical tunneling 
an outside electron is added to the Luttinger liquid, the 
plasmons spread the charge through the system, akin to 
the ripples from a raindrop on the surface of a pond. At 
zero temperature, the tunneling electron does not have 
the necessary energy to excite the plasmons. As a result, 
the tunneling conductance between a normal metal and 
a Luttinger liquid, or between two Luttinger liquids, is 
suppressed at low temperature as a power law (2,0). 

Even more interesting is the case of resonant tunneling, 
in which a single tunnel barrier between Luttinger liquids 
is replaced by a resonant level formed in a double-barrier 
quantum structure. Starting with the seminal papers by 
Kane and Fisher Q , this problem has received significant 
theoretical attention @-ll|. Perhaps the most spectac- 
ular prediction is the existence of resonance peaks with 
perfect conductance (full transparency), but vanishingly 
small width (infinite life time) at zero temperature [7|. 
These resonances require two identical Luttinger liquids 
that are symmetrically coupled to the resonant level. Al- 
though several experiments addressed resonant tunneling 
in a Luttinger liquid [12h14| |. controlling the tunneling 
strength was never attempted. 

In this work, we create a system analogous to a Lut- 
tinger liquid by properly designing electron interactions 
in the resonant level's environment and, for the first time, 
tune the system to the symmetric coupling point. In 
marked contrast with suppressed low-temperature con- 
ductance through a single tunnel barrier, we find that a 
resonant level, symmetrically coupled to two leads, re- 
tains a unitary conductance of e 2 /h, corresponding to 



perfect transparency. Simultaneously, the width of the 
resonance tends to zero as a power law of temperature. 
We argue that this unique behavior results from many- 
body physics of interacting electrons and signals the pres- 
ence of a quantum phase transition (QPT) [3, EH. QPTs 
found in strongly correlated bulk materials are often ex- 
plained by invoking interactions between local sites and 
collective modes |J|. In the same spirit, our observation 
provides the first example of a QPT in a highly tunable 
system, which emulates such a "local site" (i.e. the res- 
onant level) embedded in an interacting host. 

Both the single barrier tunneling and the double- 
barrier resonant tunneling regimes are realized in short 
(~ 300 nm) segments of carbon nanotubes (Fig. la). 
The representative electrical conductance through such a 
sample is shown in Fig. lb. The size quantization of elec- 
tron states in the nanotube, combined with the mutual 
repulsion of the electrons, result in a "Coulomb blockade" 
pattern [H,[l7|- 

We first show Luttinger liquid-like properties in tun- 
neling through a single barrier by tuning the gate voltage 
to Coulomb blockade valleys Y and Z of Fig. lb, where 
no low-energy excitations exist in the nanotube. Elec- 
trons are then transmitted through by the co-tunneling 
processes (see discussion in the Supplementary Informa- 
tion) . These processes are almost energy- independent on 
energy scales smaller than the nanotube charging energy 
and level spacing (both meV's) - the nanotube should 
behave just like a single tunnel barrier. 

The conductance in valleys Y and Z measured vs. bias 
V shows a surprising zero-bias anomaly (ZBA), which 
gets progressively deeper as the temperature decreases 
(Fig. lc). Since the shape of the ZBA in the two val- 
leys is the same up to an overall scale factor, the ex- 
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FIG. 1: Emulating Luttinger liquid with resistive en- 
vironment, a, AFM image of the sample. The nanotube 
(marked CNT) is contacted by two metal leads (S and D), 
forming a quantum dot. Two side-gates (SGi and SG2) con- 
trol the coupling of the dot to the two leads, as used later 
in Fig. 2. b, Differential conductance G = dl/dV of a sim- 
ilar sample as a function of the back-gate voltage V^atc- We 
focus on Coulomb blockade valleys Y and Z, in which elec- 
tron transport is conducted through co-tunneling processes. 
T = 1.8 K. c, Differential conductance measured vs. bias 
V demonstrates a pronounced zero-bias anomaly in both val- 
leys. Temperature: 1.7 K to 30 mK (top to bottom), d, 
Inset: zero-bias conductance G(0,T) clearly shows a power 
law dependence on T with the same exponents found in both 
valleys. Main: G(V, T) data measured in valley A at different 
temperatures can be rescaled to collapse on the same univer- 
sal curve, described by the theoretical expression of Ref. [22I ] 
(yellow line). 



istence of the ZBA is not due to the nanotube itself. 
Indeed, the distinct feature of our samples is the metal 
leads to the nanotube, which are made rather resistive 
(kfi'a). 



"Tunneling with dissipation" [18| between such 
resistive leads is known to result in suppressed conduc- 
tance dl/dV = G oc max(k B T,eV) 2r {U^M- This ex- 
pression is similar to the power-law suppression expected 
for tunneling in a Luttinger liquid; however no real Lut- 
tinger liquid is present in our sample 23j, and the ex- 
ponent r = e 2 R/h is determined by the ratio of the re- 
sistance of the leads R to the quantum resistance h/e 2 . 
Note the highly unusual appearance of the resistance in 
the exponent, which allows us to control the strength of 
tunneling suppression simply by changing R. 

Experimentally, the zero-bias conductance scales as 
G(0,T) cx T 2r , with the same exponent 2r ~ 0.6 found 
in both valleys (inset to Fig. Id); this value is consistent 
with the leads resistance (R « 6.5 kfi in this sample). 



Furthermore, we can rescale the whole set of G(V, T) 
curves measured in valley Y as shown in Fig. Id, which 
presents G(V,T)/G(0,T) as a function of eV/k B T - a 
dimensionless ratio of bias to temperature. The yellow 
curve overlaying the symbols is the result of the full the- 
oretical expression describing tunneling with dissipation 
22], in which we use the same value of r — 0.3 extracted 
from the temperature dependence. 

The expression used to fit the data in Fig. Id is iden- 
tical to the one describing tunneling between two Lut- 
tinger liquids 



24] 



The similarity may be understood 
qualitatively: both for tunneling in a dissipative environ- 
ment and in a Luttinger liquid, the tunneling electron's 
charge couples to a continuum of bosonic modes (plas- 
mons); at zero energy (temperature or bias) the electron 
cannot excite the modes, and tunneling is suppressed. 
Furthermore, the formal mapping of the two problems 
has been demonstrated for the single-barrier case in Ref. 
25j. The recipe is to replace the Luttinger interaction 
parameter g by l/(r + 1): e.g. the case of vanishing 
dissipation, r = 0, corresponds to the non-interacting 
Luttinger liquid, g = 1. It is important to realize that 
for r 7^ electrons do in fact interact with each other 
through their coupling to the bosonic modes. We use the 
analogy between tunneling in a dissipative environment 
and tunneling in a Luttinger liquid through the rest of 
this text [23j]. 

Having established Luttinger liquid-like behavior in a 
single barrier tunneling, we now turn to the main fo- 
cus of this paper: resonant tunneling between interacting 
leads. We study single-electron conductance peaks, simi- 
lar to those shown in Fig. lb, but measured on a different 
sample. A key feature of our experiment is the use of ad- 
ditional side gates to tune the coupling of the resonant 
level to the leads (Fig. la). Fig. 2a shows the differen- 
tial conductance map as a function of the side and back 
gate voltages. Clearly, the heights of the peaks change 
along the traces. For several of the peaks, the conduc- 
tance reaches a maximum at some intermediate value of 
the side gate voltage, indicating that the tunneling rates 
from the resonant level to the source and the drain are 
equal: T$ = ("symmetric coupling"). 

We focus on peak X of Fig. 2a, with the side gate tuned 
so that the tunneling is either symmetric (Fig. 2b), or 
rather asymmetric (Fig. 2c). Clearly, the two cases be- 
have in markedly different ways: in the asymmetric case, 
the peak height decreases at low temperatures, while the 
width saturates [3]. In the symmetric case, the peak 
width decreases, while the peak height grows and reaches 
e 2 /h [27l| . It is remarkable that the resonant tunneling 
conductance can reach the unitary limit despite coupling 
to the interacting leads, which suppress tunneling in the 
single barrier case. 

To account for the observed behavior, in the supple- 
mentary material we present a model of a resonant level 
connected to two electron reservoirs, with excitation of 
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FIG. 2: Resonant lineshape: symmetric and asymmet- 
ric cases, a, Zero-bias differential conductance as a function 
of V^ a te and the voltage Vsg applied to one of the side gates. 
Several peaks reach a maximal conductance of e 2 /h along 
their traces in the range shown here. The base temperature 
is T = 50 mK; perpendicular magnetic field of 6 T is applied 
to select a single spin species, b, c, Resonant conductance 
measured on the peak marked 'X' in panel a at several tem- 
peratures as a function of AV^ate, gate voltage relative to 
the center of the peak, (b) symmetric and (c) asymmetric 
coupling cases (the side gate voltages are fixed at the values 
indicated by white lines in Fig. 2a). As the temperature is 
reduced, in the symmetric case the peak becomes taller and 
narrower; in contrast, in the asymmetric case the peak be- 
comes shorter and its width saturates. 



environmental modes represented by a dynamical phase 
associated with the tunneling matrix element We 
show that the analogy between tunneling with dissipation 
and tunneling in a Luttinger liquid @, H^, [H, [29[ further 
extends to our case of resonant tunneling. Based on our 
mapping and the Luttinger liquid predictions [7-ll|, in 
the case of symmetric coupling we expect the peak height 
to saturate at e 2 /h (spinless case), and the resonance 
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FIG. 3: Resonant peak parameters at different degrees 
of asymmetry. Conductance peak height (a) and width (b) 
measured at several values of Vsg, which controls the degree 
of the tunnel barrier asymmetry. (Different peak, but the 
same sample as in Figure 2.) Note that in the symmetric 
case, the peak height saturates at e 2 /h, while the peak width 
monotonically decreases with lowering temperature. In the 
asymmetric case, the behavior is the opposite: the width of 
the peak saturates, while the peak height decreases. 



width to scale at low temperatures oc T' A r+1 ) For 
asymmetric coupling, the resonance width is predicted 
to saturate at low enough temperature 0, [l(| , while the 
peak height should scale to zero as G oc T 2r , featuring 
the same exponent as in the single barrier (non-resonant) 
tunneling. 

Our experiment clearly corroborates these predictions 
(Fig. 3). Quantitatively, we extract r ps 0.75 from the 
scaling of the asymmetric peak height, which agrees with 
the leads resistance in this sample. The width of the 
symmetric peak scales with an exponent of 0.45, consis- 
tent with r/(r + 1) ~ 0.43. (We discuss the accuracy of 
extracting the exponents in the supplementary.) Overall, 
application of Refs. f§-fll1| to our experiment describes 
the observed behavior remarkably well, both in the sym- 
metric and asymmetric cases. 

Note that the width of the conductance peak in the 
symmetric case monotonically decreases with lowering 
temperature. In the limit of zero temperature, we expect 
that the conductance will be equal to zero everywhere, 
except for a singular point at the center of the peak. 
When tunneling asymmetry is introduced, the singular 
point disappears, and the low-temperature conductance 
tends to zero at any V^ a te- This behavior indicates a 



quantum phase transition 15| for symmetric coupling, 
T S = T D . 

Indeed, in the supplementary material, we map our 
model in the r = 1 case, following Ref. 11] . onto the ex- 
otic two-channel Kondo model 30, 31 1 , for which a QPT 
is known to occur exactly for symmetric coupling In 
both models, the origin of the quantum critical behavior 
is the competition between the two channels attempting 
to screen the local site (spin or resonant level). 

The intermediate values of 1 > r > do not allow for a 
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FIG. 4: Phase diagram and the quantum critical point. 

a, Conductance in the symmetric coupling case, plotted as a 
function of temperature at different values of gate voltage, 
starting with AVgate = (top curve). Note the similarity to 
Figure 3a. b, Proposed phase diagram: the quantum critical 
point at the center (symmetric coupling and AV ga tc = 0) has 
unitary conductance. Any deviations from this point result 
in vanishing conductance at T = 0. 



at Ts — T^), AV^ a tc = [lij]. The four quadrants rep- 
resent the states of the nanotube filled with N or N + 1 
electrons, or coupled more strongly to either the source 
or the drain. The boundaries between the quadrants are 
smeared, and at T = the conductance tends to zero 
everywhere, except for the quantum critical point. 

In conclusion, we have investigated resonant tunneling 
between interacting leads emulating Luttingcr liquids. 
For symmetric coupling of the spinless resonant level to 
the two leads, and on-resonance, the low-temperature 
conductance saturates at the unitary value of e 2 /h. We 
associate this behavior with a quantum critical point, 
which exists at Ts = Tr> in the presence of a finite in- 
teraction strength r > 0. Moving away from this point 
by inducing tunneling asymmetry results in suppression 
of conductance at low temperature and smearing of the 
QPT. Our work is the first example of a QPT in a highly 
tunable system, in which many parameters can be con- 
trolled, including the strength of interactions. 



simple interpretation in terms of any Kondo model with 
non-interacting leads, but represent a continuous evolu- 
tion between the non-interacting resonant level at r = 
and the two-channel Kondo model [H, |33j • The critical 
exponents describing the system parameters close to the 
quantum critical point are not fixed, but are controlled 
by the value of r. Thus our system not only provides new 
insights into the two-channel Kondo model - a paradig- 
matic example of quantum criticality - but also gives ac- 
cess to a new family of quantum critical points for r > 0. 

The QPT observed here is different from the various 
QPTs observed " 



34] and predicted [28, 35] in quantum 
dots coupled to a single screening channel - indeed, there 
the QPTs are of the Kosterlitz-Thouless-type, while in 
our case the QPT is of second order (see the Supple- 
mentary Information). Furthermore, in our case, the key 
ingredient that enables the QPT is the symmetric cou- 
pling to the two leads, which allows for their competition; 
the interaction in the leads (finite r) prevents their hy- 
bridization. 



Finally, the conductance in the symmetric case can 
be plotted as a function of temperature for several val- 
ues of AV^atc (Fig. 4a). The similarity with Fig. 3a is 
striking: apparently, one can tune away from the uni- 
tary resonance either by inducing asymmetry (Fig. 3a) 
or by applying the gate voltage (Fig. 4a) 0] with virtu- 
ally the same results. Note that the downturn of peak 
height in Figures 3a and 4a occurs at increasingly lower 
temperature as either the degree of asymmetry or AVg atc 
is reduced. Clearly, a new energy scale is emerging in the 
system, controlled by proximity to the quantum-critical 
point 0, [t| [ljj. We anticipate that this scale should 
vanish exactly at that point. We therefore propose the 
phase diagram of Fig. 4b, with a quantum critical point 
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NANOTUBE QUANTUM DOTS 

The nanotubes are grown by chemical vapor deposi- 
tion from a CH4 feedstock gas |l| on a Si/SiOo substrate 
coated with Fe/Mo catalyst nanoparticles usually 
producing nanotubes with diameters of about 2 nm. Indi- 
vidual nanotubes are contacted by metallic leads, thereby 
forming quantum dots, and controlled by three gates: a 
back gate that changes the number of electrons in the 
nanotube and two side gates (SGi and SG2), located 
closer to either the source or the drain electrodes. Ap- 
plying the side gate voltage Vsg modifies the relative 
strength of tunneling from the nanotube to the source 
and the drain. It turns out that it is sufficient to bias 
only one of the side gates, as done in this paper. 

Figure IB of the main text shows a typical plot of the 
nanotube conductance vs. gate voltage (the "Coulomb 
blockade pattern" ) . A group of 4 peaks of similar height 
is visible on the left, corresponding to filling a 4-electron 
"shell" 0, [B[ . The peaks in the group are separated from 
neighboring groups by wider Coulomb blockade valleys Y 
and Z, in which an integer number of shells is completely 
filled. In this regime, the electron transport through the 
nanotube is determined by the co-tunneling processes 
(Figure SI), which are almost energy-independent on 
the energy scales smaller than the charging energy or 




FIG. SI: Schematic: an energy diagram of a Coulomb block- 
ade valley, in which an integer number of 4-electron shells 
is filled (like valleys Y or Z in Figure 1 of the main text). 
Cotunneling is the dominant process that contributes to the 
electron transport through the quantum dot: an electron from 
the source can virtually occupy a high energy orbital in the 
dot and then tunnel out; alternatively an electron can tunnel 
out of a filled orbital, followed by an electron from the lead 
filling the empty state. 



the level spacing (i.e. meV's). In this case, the nan- 
otube essentially behaves as a lumped tunnel junction. 
As discussed below and in the main text, this allows us 
to in situ characterize the suppression of tunneling con- 
ductance due to the resistive leads (Figures 1C-D of the 
main text). 



DISSIPATIVE ENVIRONMENT 



The strength of dissipation in the resistive leads is 



characterized by 



'R/h, where R is the total re- 



sistance of the two leads. The leads are made from a 
Cr/Au (lOnm/lnm ) film with a resistivity of 75 fi/D, 
and their total room temperature resistance is estimated 
at ~ 6.5 ki7 for the sample shown in Figure 1 of the main 
text, and ~ 17 kO for the sample shown in Figures 2-3. 
These numbers yield r ~ 0.25 for the first sample, and 
r ~ 0.65 for the second sample. However the film resistiv- 
ity increases by at least 10% at low temperature, making 
these estimates consistent with r = 0.3 we extract in Fig- 
ure ID for the first sample, and r = 0.75 extracted for 
the second sample (see the next paragraph). 

According to the theory of tunneling with dissipation 
(also referred to as "environmental Coulomb blockade" 



12j), in order to determine the dissipation strength, 



one has to consider the impedance of the whole circuit 
at frequencies corresponding to temperature or bias (i.e. 
hu> ~ fcflT or eV). In our case, this frequency range 
extends from - 1 GHz up to - 100 GHz. The litho- 
graphically made resistive leads connect the nanotubes 
to much larger pads (hundreds of microns on a side), 
whose capacitance short-circuits the high frequency fluc- 
tuations. Therefore, only the on-chip resistance of the 
leads contributes to dissipation. We also estimate that 
the distributed capacitance of the resistive leads can be 
neglected in this frequency range, so that the leads be- 
have as simple frequency - independent resistors. 

We choose to demonstrate two different aspects of the 
observed behavior in two samples. The sample shown in 
Figure 1 has the smaller dissipation strength of r = 0.3, 
resulting in the striking cusp G cx \V\ 2r in Figure 1C. In 
Figures 2-4, we extract the peak width proportional to 
T r /( r + 1 ). The relatively large dissipation strength of this 
sample is required to distinguish this non-trivial expo- 



7 



nent a from the dependence T r predicted by the lower- 
order theoretical considerations. 

Tunneling current I(V, T) through a single tunneling 
barrier in the case of Ohmic dissipation was obtained in 
Ref. [Hj]: 



I(V, T) oc VT 2 



r[r + i 



; ev 1 

J 2irk B Ti 



r[i 



eV 



'2irk B Ti 



(1) 



We numerically differentiate this expression with respect 
to V to fit the data in Figure Id of the main text. 



ESTIMATES OF THE EXPONENTS 

Asymmetric peak height. Figure 3A shows the temper- 
ature dependence of the peak height for various degrees of 
asymmetry between the tunneling rate from the level to 
the source and the drain electrodes, T$ and Td- Depend- 
ing on the asymmetry, the peak conductance exhibits a 
range of cross-over behaviors. For peaks with low degree 
of asymmetry (such that |r s -r_D| ~ k B T < |r s + r |, 
top curves), the proper temperature scaling of conduc- 
tance is not yet fully developed, so care must be taken to 
avoid extracting an incorrect exponent. Hence we study 
the conductance scaling of the most asymmetric peak of 
Figure 3A (reproduced here in Figure S2), which demon- 
strates a fully developed power-law behavior in the ac- 
cessible temperature range. A linear curve is fitted to 
the data on a log-log scale in the temperature range of 
0.08 K < T < 1.2 K. The extracted slope of 1.47 ± 0.04 
corresponds to the value of r = 0.74 ± 0.02. 
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FIG. S2: Conductance of an asymmetrically coupled peak vs. 
temperature (same data as the lowest curve in Figure 3A), 
and the power-law fit. The uncertainty in the data points is 
about the symbol size. 



Symmetric resonance width. The full-width at half- 
maximum (FWHM) of the symmetric peak (Figure 2A) 
is plotted in Figure S3. The fit shown in red yields an 
exponent of 0.43 ± 0.01, surprisingly close to the ex- 
pected value of r/(r + 1) s» 0.43. (In order to avoid over- 
estimating the exponent, the T > 1 K data points are 
not included in the fit. This higher temperature range 
corresponds to the transition to the sequential tunneling 
regime k B T > T, in which case the peak width becomes 
proportional to T .) 



> 
E 




T(K) 

FIG. S3: Full-width at half-maximum extracted from the 
symmetric peak in Figure 2B. In order to exclude the contri- 
bution from the sequential tunneling regime ksT > V, only 
the data points measured at T < 1 K are used for the fit shown 
in red (see text). The vertical axis has to be multiplied by 
the experimentally determined "gate efficiency factor" of w 
0.2, which converts AV ga te units to the actual energy of the 
resonant level (see e.g. Ref [l3|]h 

Before proceeding to the theoretical model, we note 
that our observations cannot be explained by the conven- 
tional Lorentzian expression for the resonance conduc- 
tance between non- interacting leads 14{. First, the line 



shape of the resonances is not Lorentzian. Second, even if 
one tries to describe the symmetric case by a Lorentzian 
with temperature-dependent tunneling rates Ts,d{T), in 
the asymmetric case the same expression would also yield 
a peak with a temperature- independent height and a van- 
ishing width, in marked contrast with our measurements 
in Fig. 2c. 



MODEL 

We model our experimental situation by a spinless level 
coupled to two conducting leads in the presence of an 
Ohmic dissipative environment (Figure S4). As we show 
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Vc 




FIG. S4: Schematic: a spinless quantum dot is coupled to two 
conducting leads with tunneling amplitudes Vs and Vd- The 
dot-leads system is attached to a voltage source V via R, the 
sum of the lead resistances. 



in the following, the electromagnetic modes of the en- 
vironment mediate interactions between the electrons in 
the leads, resulting in a Luttinger liquid-like behavior. 
We start with a system Hamiltonian given by 



H = H) 



Dot 



Leads 



H T + H 



Env 



(2) 



and describe the individual terms in this section, -ffoot = 
e^d'd is the Hamiltonian of the dot with the en- 
ergy level €d and the electron creation operator d) . 
#Loads = E^s^E^^L*" represents the electrons 
in the source (S) and drain (D) leads. 

Ht describes the tunneling between the dot and the 
leads: 



H T = V S I>Le-^d+) + V D 5>£ 

k 



k 



^kD^ 



'd + h.C.), 



(3) 

where the operators Ps/d represent the phase fluctua- 
tions of the tunneling amplitude between the dot and the 
S/D lead. These phase operators are canonically conju- 
gate to the operators Qs/d corresponding to charge fluc- 
tuations on the S/D junctions. This is a standard way to 
treat macroscopic quantum tunneling in the presence of 
a dissipative environment [8J] , which is valid for electrons 
propagating much slower than the electromagnetic field 

It is useful to transform to variables related to the total 
charge on the dot. To that end, we introduce Q two new 
phase operators, tp and ip, related to the phases Ps/d by 



PD = K D ip - Ip , 



(4) 



where k s / d = Cjj/g/(Cs + Co) in terms of the capci- 
tances of the two dots, C$/d- ^ is the variable conjugate 
to the fluctuations of charge on the dot Q c = Qs ~ Qd 
and so couples to voltage fluctuations on the gate which 
controls the energy of the dot's level. Likewise, ip is the 
variable conjugate to Q = (CsQd + CdQs)/{Cd + Cs)- 



Assuming Cs — Cd, we have ps — (p/2 + "4> an d 
Pd = <^/2 - ip. 

The gate voltage fluctuations can be disregarded in our 
experiment because the capacitance of the gate is negli- 
gible, C g <C C s / D . (The opposite limit of a noisy gate 
coupled to a resonant level was considered in Ref.[l6|.) In 
fact, for our purposes, the coupling of the fluctuations of 
the total charge on the dot to the environment can be 
neglected. Thus, only the relative phase difference be- 
tween the two leads remains [1, H3] , and the tunneling 
Hamiltonian becomes 



V S ^ C 



kS e 



%d + h.c.) + V D ^2{c\ D e l ^d + \\.c.) 



(5) 

The last part of Eq. Q is the Hamiltonian of the envi- 
ronment, i?Env |8l. Il8lll9||. The environmental modes are 
represented by harmonic oscillators described by induc- 
tances and capcitances such that their frequencies are 
given by uj n — l/y/L n C n . These oscillators are then 
bilinearly coupled to the phase operator ip through the 
oscillator phase: 



91 

2C 



N 

E 



Tn 

2C n 



1 

2L, 



(6) 



MAPPING TO LUTTINGER LIQUIDS 

In this section, we demonstrate the mapping of our 
model to that of a resonant level contacted by two Lut- 
tinger liquids. In carrying out this mapping, we follow 
closely previous work on tunneling through a single bar- 
rier with an environment 2(J 2ll and the Kondo effect 




in the presence of resistive leads [Lj 

The two metallic leads in our case can be reduced to 
two semi-infinite one-dimensional free fermionic baths, 
which are non-chiral 12211 . By unfolding them, one can 
obtain two chiral fields [22| , which both couple to the dot 
at x — 0. We bosonize the fermionic fields in the standard 

-^=F S/D exp[i(p s/D (x)]. Here, 
Ps/d is the bosonic field, and a 



mm c s/D (x) 



way 

Fs/d is the Klein factor, 



is the short time cutoff. Defining the flavor field <p>f and 
charge field <p c by 



4> } 



v/2 ' YC ~ V2 
we rewrite the Hamiltonian of the leads as 
v F 



H 



Leads 



4-7T 



d.r 



{d x <pcf + [{d x (p f f 



(7) 



(8) 



The tunneling Hamiltonian then becomes 
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Ht — Vs 
+ V D 



1 



1 

/2ita 



F$ exp 
Fjj exp 



.Mo) + 4>f(P) 



V2 . 

.Mo) -Mo) 



V2 



e l ~d + h.c 
e l %d + h.c 



(9) 



Note a key feature of Ht'- the fields <p and </>/(0) enter in the same way. Thus we wish to combine these two fields, a 
process which will lead to effectively interacting leads as in a Luttinger liquid. 

To carry out such a combination, since the tunneling only acts at x = 0, it is convenient to perform a partial trace 
in the partition function and integrate out fluctuations in <t>ri f (%) f° r au x away from x = [24| . For an Ohmic 
environment, one can also integrate out the harmonic modes 17, la, H3- Then, the effective action for the leads and 
the environment becomes 



Qcff 

°Lcads+Env 



1 / R 

p J2 l w »l (l^K)l 2 + I0/K)l 2 + 



(10) 



where Rq — h/e 2 , R is the total resistance of the leads, and Lu n = 2irn/ '(3 is the Matsubara frequency. In this discrete 
representation, it is straightforward to combine the phase operator (p and the flavor field <pf, in order to maintain 
canonical commutation relations while doing so, we use the transformation 



7 = \f9f [ <f>f + -j^f 



R Q Pf V R ^ 



(11) 



where gf = 1/(1 + R/Rq) < 1. Now, the effective action for the leads and environment (excluding tunneling) becomes 

(12) 



ocff 

^Lcads+Env 



= ~ E (I^K)I 2 + l^/K)l 2 + I^K)| 2 ) , 



while the Lagrangian for the tunneling part reads 



Lit — —Vs 



F S 
f2na 



-Vn 



D 



>2na 



(13) 



(Here, we have formally introduced the parameter g c = 1 
to describe the noninteracting field <f> c .) Indeed, we see 
that the phase (p has been absorbed into the new flavor 
field (f>'f at the expense of a modified interaction parame- 
ter gf, while the new phase fluctuation ip' decouples from 
the system. In what follows, we will drop the prime from 
the operator <f>'f. 

It turns out that one obtains a very similar effective 
action by starting from a model with resonant tunneling 
between Luttinger liquids 13, HI- The two models are 
equivalent if the interaction parameters g c and gf in our 
model are made equal to the single interaction parameter 
g of Ref. l2J. Similar mappings were obtained for a spin- 
ful model in the absence of charge fluctuations (Kondo 
regime) in Ref. 17 and for a dissipative dot coupled to a 



single chiral Luttinger liquid in Ref. [21 



As discussed in the main text, the r = case (gf = 1) 
reduces to a resonant level without dissipation, while the 
r = 1 case (gf — 1/2) can be mapped onto the two- 
channel Kondo model. To understand the relation be- 
tween our model and the two-channel Kondo model, we 



apply a unitary transformation 26, 27]], U = exp[i(cft d — 
l/2)0 c (O)/v / 2], to eliminate the <p c field in the tunneling 
Lagrangian, Eq. (|13l) . At the same time, an extra elec- 
trostatic Coulomb interaction between the leads and the 
dot is generated. Introducing in addition a bare electro- 
static Coulomb interaction ([/&), we have 



Hr 



V2 



(U b - l)(dU - l/2)d x Mx = 0) . (14) 



For the special value gf — 1/2, we can refermionize the 
problem by defining ip c j = e l ^ c - f / \j2na. Then, at Vs = 
Vd (and r — 1), our model is mapped onto the two- 
channel Kondo model 26-2^], which shows non- Fermi- 
liquid behavior [30]. At the Toulouse point (Ub = 1 so 
that He = 0), the model reduces to a Majorana resonant 



level model [26j, |30|. Finally, for r close to 1 ( 



our case), one can similarly map our system onto a two- 
channel Kondo model with (interacting) Luttinger liquid 
leads 2|| 2t|. The effective interaction parameter g a in 
this case is determined by the residual 1 — r. 



10 



QUANTUM PHASE TRANSITION AND 
SCALING RELATIONS 

Having established the relation between our problem 
and Luttinger liquid physics, we can draw on the very 
extensive theoretical work concer ning re sonant tunnel- 
ing in a Luttinger liquid 2J, |27|, [31-35] to reach con- 
clusions about the scaling and phase transitions implied 



by our model. To understand how the interacting envi- 
ronment affects the low temperature physics, it is conve- 
nient to rewrite the model, following Refs. 24 andl36l in 
the "Coulomb-gas" representation, which can be accom- 
plished by expanding the partition function in powers of 
Vs and Vd- After integrating out 4>f( T ) an d <I>c(t) in each 
term, one obtains a classical one-dimensional statistical 
mechanics problem with the partition function 



z = ^ ^2 y g Ei( 1 +9-P-)/2 v -E i (l-9 i Pi)/2 

a=± n { Qi= ±} 

x J dr 2n J dT 2n - 1 J dr\ exp | ^ V tJ } 

*<» 

l [n — T- 

Vij = - — [qiqj + K 1 p l p j + K 2 (p l q 3 + Pjqi)} In 

2 9f \ t c 



exp < ed 



.1 - o \ - 

p— — 1~ & 2^ Pin 

KK2n 



(15) 



We consider the on-resonance case, ed = 0, so that the 
last term in the partition function is equal to 1. There are 
two types of charges in this ID problem 24]: qi charge 
and pi charge, both of which can be ±1. The total system 
is charge neutral, J^. qi — ^2iPi = 0. Physically, the qi 
charge corresponds to a tunneling event between the dot 
and the leads (+1 : to the right (D); -1: to the left (S) 
). The pi charge corresponds to hopping onto (+1) or off 
(-1) the dot. 

Following the method developed for resonant tunnel- 
ing between Luttinger liquids in Ref. Q, one obtains the 
renormalization group (RG) equations and the phase di- 
agram. The sign of the pi charge must alternate in time, 
while the qi charge can have any ordering satisfying the 
charge neutrality constraint. Therefore, the interaction 
between the qi charges does not vary in the RG flow, 
while the interaction strength between pi charges, K\, 
does renormalize. The bare K\ (initial value in RG flow) 
is K^ arc = gf/g c ', in fact, by taking into account the 
coupling of the fluctuations of the total dot charge to 
the environment — an effect neglected here — one can show 
that K\ &ve = 1 [37]. The bare value of K 2 is zero, but it 



(1-0) 



(1,1) 



will be generated in the RG flow. These are the same con- 
ditions as in resonant tunneling between two Luttinger 
liquids . Hence, our model is mapped onto the model 
of resonant tunneling between Luttinger liquids, with the 
RG equations for the interaction strengths and tunneling 
couplings as in Ref. [24j. 

The resulting schematic RG flow diagram in the Vs-Vd 
plane, valid for r < 1 and on-resonance, is shown in Fig- 
ure S5. A remarkable feature of this RG flow diagram is 
that a second order quantum phase transition can be re- 
alized by tuning the tunneling matrix element across the 




(0,1) 



FIG. S5: Schematic representation of the renormalization 
group flow for the two tunneling amplitudes, Vs and Vd [H[ • 
The diagonal corresponds to the symmetric coupling case, and 
flows into the strongly coupled quantum critical point at (1,1) 
which corresponds to a single homogeneous Luttinger liquid. 
Point (1,0) describes the strong coupling point for Vs, while 
Vd = 0; similarly, point (0,1) describes the strong coupling 
point for Vd, while Vs = 0. 



symmetric coupling line, Vs = Vd. Indeed, this critical 
value separates the flows terminating in the stable fixed 
points denoted as (0, 1) and (1, 0) in Figure S5. These two 
phases correspond to the dot merging with the D lead, 
while the S lead decouples, or vice versa. The unstable 
fixed point denoted as (1,1) is reached starting exactly 
at Vs = Vd', it corresponds to a homogeneous Luttinger 
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liquid which therefore has conductance G = e 2 /h [24, 31 1 . 

The proximity to the strongly coupled point (1,1) de- 
termines the quantum critical behavior and the critical 
exponents observed in the experiment. By analyzing the 



universal scaling function [24j , the width of the resonant 
peak is found to scale to zero as V oc T 1 ^ g f = T' r /( 1 + r ). 
For the case of asymmetric tunneling Vs ^ Vb, either 
Vs or Vd flows to zero so that the height of the resonant 
peak scales as G oc r 2 ( 1 /s/- 1 )= ^_T 2r El, while the 

" 11 S 



27 



width of the resonance saturates 

These power-law dependencies are typical for the crit- 
ical behavior near a second-order phase transition. We 
would like to stress that this behavior qualitatively dif- 
fers from the Kosterlitz-Thouless (KT) type of transition, 
commonly encountered in quantum impurity models [38[ . 
For example, several recent theoretical works predict KT 
transitions in quan tum dots coupled to a single inter- 
active lead 39|-43j]. There, the QPT occurs when the 
dissipation exceeds a certain critical value. In our case, 
the crucial ingredient that enables the QPT is the sym- 
metric coupling to the two leads, which allows for their 
competition, while the dissipation strength can be rela- 
tively weak (but greater than zero). 
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